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Abstract 



After having investigated the regular prisms and prism tilings in SL2R 
space in previous work [ 13] of the second author, we consider the problem 
of geodesic ball packings related to those tilings and their symmetry groups 
pq2i. SL2R geometry is one among the eight Thurston geometries that 
;J] , can be derived from the 3-dimensional Lie group of all 2 x 2 real matrices 

. .' with determinant one. ^-^_^ 

In this paper we consider geodesic spheres and balls in SL2R space 
(even in SL2R), if its radius p € [0, |) and determine their volumes. More- 
over, we consider the prisms of the above space and compute their volumes, 
define the notion of the geodesic ball packing and its density. We develop 
a procedure to determine the densities of the densest geodesic ball packings 
for the regular prism tilings, moreover for their group pq2i, and apply this 



*Mathematics Subject Classification 2010: 52Cn^52C22, 52B15, 53A35, 51M20. 
Key words and phrases: Thurston geometries, SL2R geometry, ball packings, tiling, regular 
prism, volume in SL2R space. 
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algorithm to them. Surprisingly we found a very dense packing in SL2R 
(Theorem 4.1), denser than those in the hyperbolic space H'^. 

We use for the computations and visualization of SL2R space its pro- 
jective model introduced by the first author in |,4J. 



1 On SL2R geometry 



The real 2x2 matrices f 1 with unit determinant ad — bc = 1 constitute a Lie 

\c a) 

transformation group by the usual product operation, taken to act on row matrices 
as on point coordinates on the right as follows 

{z\ z^) ("^ ^\ = {z^d + z'c, z% + z'a) = («;°, w'). (1.1) 

This group is a 3-dimensional manifold, because of its 3 independent real coordi- 
nates and with its usual neighbourhood topology [|8l, [fTSl . In order to model the 
above structure in the projective sphere VS^ and in projective space V^ (see [|41) 
we introduce the new projective coordinates (x°, x^, x"^, x^) where 

CX • — •Jj I iX> . C/ . — tiy I iX^ . C- • — tiy I •Aj . tX . — tiy Uj « 

with positive, then the non-zero multiplicative equivalence as a projective free- 
dom in VS^ and in V^, respectively. Meanwhile we turn to the proportionality 
SL2R~PSL2R, natural in this context. Then it follows that 

0>bc~ad= -a;°x° - x'^x^ + x^x^ + x^x^ (1.2) 

describes the interior of the above one-sheeted hyperboloid solid 'K in the usual 
Euclidean coordinate simplex, with the origin i?o(l; 0; 0; 0) and the ideal points of 
the axes E'^{^; 1; 0; 0), S|°(0; 0; 1; 0), E|°(0; 0; 0; 1). We consider the coUineation 
group G* that acts on the projective space V^ and preserves a polarity, i.e. a 
scalar product of signature ( 1-+), this group leaves the one sheeted hyper- 
boloid solid "H invariant. We have to choose an appropriate subgroup G of G^, as 
isometry group, then the universal covering group and space T-LofH will be the 

hyperboloid model of SL2R (see (4]). 
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The specific isometries S(0) (0 G R) constitute a one parameter group given 
by the matrices 



S(0) : (sM) 



( COS0 sin0 \ 

— sin cos 

cos — sin <; 

y sin0 COS0 j 



(1.3) 



The elements of S(0) are the so-called ^Z?re translations. We obtain a unique fibre 
line to each X(x°; x^; x^; x^) G "H as the orbit by right action of S(0) on X. The 
coordinates of points lying on the fibre line through X can be expressed as the 
images of X by S (</>): 



/ 1 2 3n ^\'i'\ ( 



x^ sin (/); x^ sin + x^ cos < 



x^ cos f^ + x^ sin 6: —x^ sin r/) + x^ cos ( 



(1.4) 



The points of a fibre line through X by usual inhomogeneous Euclidean coordi- 
nates X = 4, y = 4, 2; = 2^, x° 7^ are given by 



(l;x;y;z) 



1; 



X + tan (f) y + z tan -z — y tan 1 



xtani 



1 — xtan0 1 — xtan( 



(1.5) 



for the projective space V^, where ideal points (at infinity) conventionally occur. 
The TV periodicity of the above maps can be seen from the formula (1.5) e.g. if 
— I < < I then —00 < X < 00. In (1.3) and (1.4) we can see the 2% periodicity 




Figure 1 : The hyperboloid model 
of (j), moreover the (logical) extension to G R, as real parameter, to have the 
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universal covers H and SL2R, respectively, through the projective sphere VS^. 
The elements of the isometry group of SL2R (and so by the above extension the 
isometries of SL2R) can be described by the matrix (a]) (see [SI and 15]) 






\±al 
-ia'of - (air + {air + (alf 



±< 



Ta', 



±al 



4\ 



Tal 



where 



Tal ±al) 



-1, 



,0\2 



0^1 



a\2 



+ {al? + (a: 



3\2 



,1^0 



— aQa2 — ^0^2 + '^o'^2 + '^0*^2 = = — aQa2 + clqC^ 



alal + alal. 



(1.6) 

We define the translation group Gt, as a subgroup of the isometry group of 
SL2R, the isometries acting transitively on the points of "H and by the above ex- 
tension on the points of "H. Gt maps the origin £"0(1; 0; 0; 0) ontoX(x°; x^; x^; x^). 
These isometries and their inverses (up to a positive determinant factor) can be 
given by the following (tl) and T^ = {tf)^^ matrices: 



T : (ti) 



(T' 



( 



X 



V 



X 

x^ 



f x^ 
x^ 



X 



X 



X 

x^ 



X 



X 



\ 



X 



\—x^ 



-X 



X 



2 '^ 

/yj-" ry'--' 



X 



-X 



-X 

-x^ 





x 
x^ 



X 



J 

—x^\ 

I 



X 
—X 



(1.7) 



The rotation about the fibre line through the origin _Eo (1; 0; 0; 0) by angle a; (— vr < 
uj < Tl) can be expressed by the following matrix (see [4]) 



ReoM: {rl{E,,uj)) 



/I \ 

10 

COS a; sin a; 

yo — sin UJ cos ujj 



(1.8) 



and the rotation Rx(w) with matrix : {rl{X,uj)) about the fibre line through 
X{x^; x^; x^; x'^) by angle cu can be derived by formulas (1.7) and (1.8) where 
Rx(w) = T^^IlEoi^)'^- Thus the above rotation Rx((^), with a specific X 
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(cosh r, 0, sinh r, 0) ~ (1, 0, tanh r, 0) has the important matrix 

( l+sinhV- „ -sinh2r- ^ 

2 sinh rsinw -. — ^ sinh 2r sin w 

— sinh r cos to -1,0 

smn 2r cos co 

' 1 

2 1 + sinhV- , --sinh2r+ 

— sinh r sin to ^ ~ 2 ^^^^ ^^ ^^^ '^ 1 
-sinh rcosw + i- sinh 2r cos ^ 

1 ^ 

--sinh2r-+ i_cosh2r+ 2 

-. — 2 sinh 2r sin w 2 cosh rsinw 

+-sinh2rcosa; +^°^^ ^^"^'^ 

1 , 

-smh2r- i_cosh2r+ 

— 2 sinh 2r sin cj . —cosh rsinw 2 

y --sinh2rcosa; +cosh rcosa; ^ 

(1.8') 
Horizontal intersection of the hyperboloid solid 'H with the plane EqE^E^ pro- 
vides the base plane of the model "H = SL2R. The fibre through X intersects the 
base plane z^ = x = in a foot point 

Z{z^ = a;°2;° + x^x^; z^ = 0; z'^ = x^x^ - x^x^] z^ = x^x^ + x^x'^). (1.9) 

We generally introduce a so-called hyperboloid parametrization by flU as follows 



x° = cosh r cos ^, 
x^ = cosh r sin (h, 



X = sinh r cos 



(1.10) 



') 



X = sinh r sin 1 



where (r, 6) are the polar coordinates of the base plane and (j) is just the fibre 
coordinate. We note that 

—XX —XX + X X + XX = — cosh r + sinh r = — 1 < 0. 

The inhomogeneous coordinates in (1.1 1), that will play an important role in the 
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later E^-visualization of the prism tilings in SL2R, are given by 

X = —r = tan0, 
x" 

x"^ , cos (9 — 6) 
^= =tanhr ^— ^, (l.H) 

x^ cos (p 

x^ , sin (6 — (j)) 

z = —r = tannr ■ . 

x^ cos (p 

The infinitesimal arc-length square can be derived by the standard method, 
called pull back into the origin. By acting of ( 1 .7) on the differentials (dx° ; dx^ ; dx"^ ; dx^ ) 
we obtain by flU, |l2l and Q that in this parametrization the infinitesimal arc- 
length-square at any point of SL2R is the following: 

{dsf = (dr)2 + coshVsinhV(ci^)2+ [{dcj)) + smh^ r{de)f . (1.12) 



Hence we get the symmetric metric tensor field gij on SL2R by components: 

/I \ 

Qij := I sinh r(sinh r + cosh r) sinh r | , (1-13) 

and 



dV = \/det{gij) dr d9 d(j) = — sinh(2r)(ir d6 d(j) 

as the volume element in hyperboloid coordinates. The geodesic curves of SL2R 
are generally defined as having locally minimal arc length between any two of 
their (near enough) points. 

By (1.13) the second order differential equation system of the SL2R geodesic 
curve is the following: 

r = sinh(2r)^0 H — (sinh(4r) — sinh(2r))^^, 

4> = 2rtanh(r)(2sinh^ {r)9 + (j)), (1.14) 

We can assume, that the starting point of a geodesic curve is the origin (1, 0, 0, 0), 
because we can transform a curve into an arbitrary starting point. Moreover, 
r(0) = 0, 0(0) = 0, 6(0) = 0, r(0) = cos(a), 0(0) = sin(a) = -^(0) in 
Table 1 for the solution of (1.14) and so the unit velocity will be achieved. 
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Table 1 


Types 




TT 

< a < - 

4 

(H^ — like direction) 




r{s,a) = arsinh ( — sinh(s-\/cos 2a) ) 

V Vcos 2a / 


fii^fi n) = — arrtan tanh^.syrns yre ) 


^ Vcos 2a ' '^ 
0(s, a) = 2 sin as + 9{s, a) 


vr 
"=4 
(light direction) 


r{s, a) = arsinh ( — s) 

9{s,a) = — arctanf-— s) 
(j){s,a) = V2s + 9{s,a) 


TT TV 

- <a< - 

4 - 2 

(fibre — like direction) 


/ cos a \ 


r{s,a) = arsinh ( """ sin(s V— cos 2a) 

Vv-cos2a / 


6(s,a) = — arctan "'" "^ — tan(sV— cos2a) 
V a/— cos 2a / 

(f){s, a) = 2 sin as + 9{s, a) 



The equation of the geodesic curve in the hyperboloid model has been de- 
termined in [|2l, with the usual geographical sphere coordiantes (A, a), as longi- 
tude and altitude, respectively, from the general starting position of (ITO), (IT 1), 
(— TT < A < vr, — I < a < |), and the arc-length parameter < s G R. The Eu- 
clidean coordinates X{s, A, a), Y{s, A, a), Z(s, A, a) of the geodesic curves can 
be determined by substituting the results of Table 1 (see dZJ) into the equations 
(ITO) and (ITl) as follows: 



Y{s, A, a) 



X{s, A, a) = tan (0(s, a)) 
tanh (r(s, a)) 



cos (0(s, a) 



cos 



[9{s, a) — (f){s, a) + A] . 



(1T5) 



Z{s, A, a) 



tanh (r(s, a)) 
sin (0(s, a) 



cos 



["( 



s, a — 0[s, a 



+ A] 
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2 Geodesic balls in SL2R space 

Definition 2.1 The distance d{Pi, P2) between the points Pi and P2 is defined by 
the arc length of the geodesic curve from Pi to P2. 

The numerical approximation of the distance d{0, P), by Table 1 and (1.15) for 
given P{X, Y, Z) from the origin O, will not be detailed here. 

Definition 2.2 The geodesic sphere of radius p (denoted by Sp^ (p)) with the cen- 
ter at point Pi is defined as the set of all points P2 in space with the condition 
d{Pi, P2) = p. Moreover, we require that the geodesic sphere is a simply con- 
nected surface without selfintersection in SL2R. 

Definition 2.3 The body of the geodesic sphere of centre Pi and of radius p in 
the SL2R space is called geodesic ball, denoted by Bp^(p), i.e., Q G Bp^{p) iff 
0<d{Pi,Q) <p. 

Remark 2.1 We choose 0(1, 0, 0, 0) as centre of the sphere and its ball by the 
homogeneity 0/SL2R. 

Fig. 2. a shows a geodesic sphere of radius p = 1.3 with centre O and Fig. 2.b 
shows its intersection with the (x, z) plane. From (1.15) it follows that S{p) is a 

simply connected surface in E^ and SL2R, respectively, if p G [0, |). Ifp > | 
then the universal cover should be discussed. Therefore, we consider the geodesic 
spheres and balls of radius p G [0, |). These will be satisfactory for our prism 
tilings. 



2.1 The volume of a geodesic ball 

The volume formula of the geodesic ball B{p) of radius p follows from the metric 
tensor gij. We obtain the connection between the hyperboloid coordinates (r, 6, 0) 
and the geographical coordinates (s, A, a) by the Table 1 and by the equation 
(1.15). Therefore, the volume of the geodesic ball of radius p can be computed by 
the following integral: 
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Theorem 2.1 



Figure 2: a, b 



Vol(B{p)) = - sinh(2r) dr dO dcp 
J B 2 



fp fi I 
= 4:11 / - smh{2r {s,a))\Ji\ da ds+ 
Jo Jo 2 

fp [^ I 
+47r / / -smh{2r{s,a))\J2\ da ds 
Jo J^ 2 



(2.1) 



where \Ji 

ds da 

corresponding Jacobians. 



dr dr 
ds da 



and similarly \ J2 \ (by Table 1 and 



d\ 



1) are the 



The complicated formulas above need numerical approximation by computer. In 
Table 2 we have collected the volumes of the geodesic balls for some radii (Fig. 3). 
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I I I I III I I I I I I I I I I I I I I I I I I I I I 

0.0 0.25 0.5 0.75 1.0 1.25 1.5 
rho 



Figure 3: The increasing function p i-)- Vol{B{p)). 



Table 2 


p Vol{B{p)) 


p = 0.1 ^ 0.004202769 


p = 0.5 ^ 0.568541467 


p = 1 ^ 5.761162741 


p = 1.3 ^ 15.54409646 


p^l ^ 34.17279403 



3 Regular prisms and prism tilings in SL2R space 

In the paper [fT3l we have defined and described the regular prisms and prism 
tilings in SL2R space. These will be summarized in this section. 

Definition 3.1 Let "P* be an infinite solid that is bounded by certain surfaces in 
the model space 0/SL2R that will be determined (later on) by ,, side fibre lines" 
passing through the vertices of a regular p-gon (V'') lying in the ,,base plane". 
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The images of solids V by SL2R isometries are called infinite regular p-sided 
SL2R prisms. Here regular means that the side surfaces are congruent to each 
other under SL2R rotations about a fiber line (e.g. through the origin). 

The common part of V^ with the base plane is the base figure of P* that is denoted 
by V and its vertices coincide with the vertices of V^ but V is not assumed to be 
a polygon. 

Definition 3.2 A bounded regular p-sided prism in SL2R is analogously defined 
if the face of the base figure V and its translated copy V* by a fibre translation by 
(1.3) and so (1.5) are also introduced. The faces V and P* are called cover faces. 

Remark 3.1 All cross-sections of a prism ,,parallel" (the intersecting planes are 
generated by SL2II fibre translations from the base plane) to the base faces are 
congruent. Prisms are named for their base, e.g. a prism with a pentagonal base 
is called a pentagonal prism. 

In [fT3l we have considered regular prism tilings 7^(1?) by prisms Vp{q) where q 
pieces regularly meet at each side edge by g-rotation. 
The following theorem have been proved: 

Theorem 3.1 There exist regular infinite prism tilings Tp{q) in SL2R/or each 
3 < p G N where -^ < g G N. For bounded prisms, these are not face-to-face. 

In this paper we consider the above regular bounded prism tilings Tp{q). Then 

each vertex of a tiling is a proper point of SL2R, and the prism is a , , topological 
polyhedron" having at each vertex one ,,p-gonal cover face" (it is not a polygon at 
all) and two skew quadrangles which lie on certain side surfaces in the model. Let 
Vp{q) be one of the tiles where V (and so V^ as well) is centered at the origin with 
vertices A1A2A3 ... Ap in the base plane. It is clear that the side curves CA^Ai+i 
(z = 1 . . .p, Ap^i = Ai) of the base figure are derived from each other by — 

rotation about the vertical x axis, so there are congruent in SL2R sense. The 
corresponding vertices B1B2B3 . . . Bp are generated by a fibre translation r given 
by (1.3) with parameter < $ G M \ {0}. The fibre lines through the vertices 
AiBi are denoted by /j, (i = 1, . . . ,p) and the fibre line through the "midpoint" 
H of the curve ca^a is denoted by /q. This /o will be a half-screw axis of ab as 
follows. 
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B=A^ 




Figure 4: 



The tiling Tp{q) is generated by its discrete isometry group Tp{q) = pq2i 

C /som(SL2R) which is given by its fundamental domain AiA20AlAlO^ a 
,, topological polyhedron" and the group presentation (see Fig. 4 for p = 3 and 

my. 



pq2i = {a, b, s : a^ = b'^ = asa ^s ^ = babs ^ 
= {a, b : aP = b-^ = ababa^b^a^b^ = 



1} 



!}• 



(3.1) 



Here a is a p-rotation about the fibre line through the origin (x axis), b is a g- 
rotation about the fibre line trough Ai and s = bab is a screw motion s : OA1A2 — )■ 
O^BpBi. All these can be obtained by formulas (1.8) and (1.8'). Then we get the 
second presentation in (3.1), i.e. abab = baba =: r is a fibre translation. Then 
ab is a 2i half-screw motion about /o = HH'^ (look at Fig. 4) that determines 
the fibre tamslation r above. This group in (3.1) surprisingly occurred in § 6 of 
our paper [7J at double links Kp^g. The coordinates of the vertices A1A2A3 . . . Ap 
of the base figure and the corresponding vertices B1B2B3 . . . Bp of the cover face 
can be computed for all given parameters by 



tanh(OAi) = b:-- 



- tan - tan - 
p 1 

1 + tan ^ tan ^ ' 
g 1 



(3.2) 
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Moreover, the equation of the curve caiA2 can be determined as the foot points 
(see (1.4) and (1.9)) of the corresponding fibre lines. For example, the data of 
Vs{q) for some 6 < g G N are collected in Table 2 by Maple computations. 



Table 2 


{p,q) 


b 


(3,7) 


^ 0.30007426 


(3,8) 


^ 0.40561640 


(3,9) 


f« 0.47611091 


(3,10) 


^ 0.50289355 


(3,50) 


^ 0.89636657 


(3, 1000) 


^ 0.99457331 


(3,oo) 


1 



3.1 The volume of the bounded regular prisms Vp{q) 

The volume formula of a sector-like 3 -dimensional domain Vol{D{^)) can be 
computed by the metric tensor gij (1.13) in hyperboloid coordinates. This defined 
by the base figure D lying in the base plane (see Fig. 5) and by fibre translation r 
given by (1.3) with parameter $ G M \ {0}. 

Theorem 3.2 Suppose we are given a sector-like region D illustrated in Fig. 5. a 
so a continuous function r = r{6) where the radius r depends upon the angle 6. 
The volume of domain D{^)) is derived by the following integral: 



Vol{D{<^))= I -smh{2r)drde 
' D 2 



* /•6»2 prie) 



Je 



sinh(2r) dr dO 



$ 



02 



(3.3) 



(cosh(2r(6')) - 1) dO. 



Let Tp{q) be a regular prism tiling above and let Vp{q) be one of its tiles. We get 
the following 

Theorem 3.3 The volume of a bounded regular prism Vp{q) [3 < p G N, -^ < 
q EN) can be computed by the following formula: 



Vol{Vp{q)) = Vol{D{p,q,^)) ■ p, 



(3.4) 
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where Vol{D{p, q, $)) is the volume of the sector-like 3 -dimensional domain that 
is given by the sector region OA1A2 C V (see Fig.S.b) and by ^ = AiBi, the 
SL2R height of the prism, depending on p, q. 





a. b. 

Figure 5: The sector region D 



4 Ball packings to regular prism tilings 

Sphere packing problems concern arrangements of non-overlapping equal spheres, 
rather balls, which fill a space. Space is the usual three-dimensional Euclidean 
space. However, ball (sphere) packing problems can be generalized to the other 
3-dimensional Thurston geometries. But sometimes a difficult problem is - simi- 
larly to the hyperbolic space - the exact definition of the packing density. In [fT4l 
we extended the problem of finding the densest geodesic ball packing for the other 
3-dimensional homogeneous geometries (Thurston geometries). In this paper we 

study the problem in SL2R geometry and develop a procedure for regular prism 
tilings and their above group pq2i in (3.1). 

We consider the above regular prism tilings Tp{q). 

The greatest ball of centre K inside of Vp{q) is denoted by B^, moreover let 
p{K) denote the radius of Bk- The images of Bk in each tile of 7^(g) form a ball 
packing Bk- 



Volumes and geodesic ball packings to the regular prism tilings... 15 

Definition 4.1 The density of ball packing Bk is 

VoI{Bk) 



5{K) 



Vol{VM) 



Our problem is to find a point K E 'Pp{q) and so the ball packing Bk where its 
density 6{K) is maximal. In this case the ball packing Bk is said to be locally 
optimal. 

Our aim is to determine the maximal radius p(Kopt) of the balls, and the max- 
imal density S{Tp{q)) = 5{Kopt)- The considered prism tiling does not have free 
parameters, thus 

5{Vp{q)) = uiM5{K)). (4.1) 

K 



4.1 The optimal ball packings to prism tilings Tp{q) 

Let Tp{q) be a regular prism tiling and let Vp{q) be one of its tiles which is given 
by its base figure V that is centered at the origin with vertices AiA2A'i . . . Ap 
in the base plane of the model (see Fig. 7) and the corresponding vertices B1B2 
B^ . . . Bp and C1C2C3 ... Cp are generated by fibre translations r and — r given 
by (1.3) with parameter $ G M \ {0}, as above. 

It is clear that the optimal ball Bk has to touch some faces of the correspond- 
ing prism Vp{q) and it is centered at the origin. The Vol{Vp{q)) is given by the 
parameters p and q and $ > p{Kopt)- The images of Vp{q) by our discrete group 

pq2i covers the SL2R space without overlap. For the density of the packing it 
is sufficient to relate the volume of the optimal ball to that of the solid Vp{q) (see 
Definition 4.1 and (4.1)). 

The maximal radius p{Kopt) of the balls can be determined for all possible 
parameters as the distance p between the origin and c^^^j (see Fig. 6), related to 
the half height $ = | — - — - of the prism, also determined by p and q (see 
and [ini). 

The above locally densest geodesic ball packings can be determined for all 
regular prism tilings Tp{q) {p, q as above). For example, in Table 3 we have sum- 
marized some results to tilings 73(g) (6 < g G N). 
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Figure 6: The optimal half prism A1A2A3B1B2B2, with optimal half sphere for 



parameters p = 3, q = 7 with the maximal radius p{Kopt) = ^ 



<p. 



Table 3 


{p,q) 


p{Kopt) 


VoI{Bk) 


Vol{VM) 


5{Kopt) 


(3,7) 


0.074800 


0.001756 


0.016785 


0.104635 


(3,8) 


0.130899 


0.009449 


0.051404 


0.183816 


(3,9) 


0.174533 


0.022497 


0.091385 


0.246177 


(3, 10) 


0.209440 


0.039048 


0.131595 


0.296731 


(3,50) 


0.272672 


0.087043 


0.636918 


0.136663 


(3, 1000) 


0.274648 


0.088981 


0.812627 


0.109498 



4.2 Multiply transitive optimal ball packings under pq2i 

We study only one case of the multiply transitive ball packings where the funda- 
mental domain of the SL2R space group pq2i is not a prism. Let the fundamental 
domains be Dirichlet- Voronoi cells where their centers are images of K = O. 
These locally densest geodesic ball packings can be determined for all possible 

fixed parameters p and g ( 3 < p G N, -^ < g G N) . The optimal radius p{Kopt) 
{K = O) is 



piK 



opt) 



mm 



artanh(Oyli), $ 



TV n 



n 



d{0, O^*^) 



p 



(4.2) 
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a 



Figure 7: The optimal cell and ball for parameters p = 3 and q = 7. 

where d{0, O^^) is the distance between O and O^^ by Definition 2.1. 

In Table 4 we have summarized some results and we noticed that in the case 
7(5(9) we have obtained a very dense ball packing whose density is larger than 
~ 0.85327609 the universal hyperbolic 3-dimensional density upper bound by 
Boroczky and Florian ^. Hence we have as a corollary 



Table 4 


ip,q) 


p{Kapt) 


VoI{Bk) 


Vol{V,{q)) 


KKopt) 


(3,11) 


0.237999 


0.057543 


0.169931 


0.338626 


(3,12) 


0.261799 


0.076892 


0.205617 


0.373960 


(3, 13) 


0.279134 


0.093489 


0.238467 


0.392044 


(3,14) 


0.287083 


0.101857 


0.268561 


0.379271 


(3,50) 


0.350810 


0.188371 


0.636918 


0.295754 


(3, 1000) 


0.370822 


0.223543 


0.812627 


0.275087 












(6,8) 


0.654498 


1.350812 


0.640122 


0.823604 


(6,9) 


0.692287 


1.624770 


1.877243 


0.858262 



Theorem 4.1 The optimal ball packing to the SL2R group pq2i above for p 
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6, q = 9 with ball centres in the p-gonal rotation axis, is denser than any ball 
packing, related to its Dirichlet-Voronoi partition, in the hyperbolic space H^. 

Remark 4.1 In a forthcoming paper we will study the above locally optimal 
geodesic ball packing arrangements and their densities for all possible param- 
eters. 

Our projective method gives us a way of investigating the Thurston geome- 
tries. This suits to study and solve similar problems, (see e.g. 161, flU, [[T2l . 

a4j). 
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